close , clear, clc;
lanbta1 = 1; h1 = 1; S1 = 1; S = 1;
Nt = 50; dt = 0.001; dh = 2;
u1 = 1; u2 = 1.001; u3 = 1; v1 = 1; v2 = 1; v3 = 1;
i = sqrt(-1);
for ii = 1:100
    T(ii) = 10 ^ ((ii - 20) / 10);
   log10_p(ii) =  2 * log10(sqrt(T(ii)) ) - 2* log10((u2+v2) / 365);
    pT(ii) = (365 * sqrt(T(ii) ) / S )^2;
end
figure(1)
plot(pT)
figure(2)
plot(log10_p)


a = ['asdf','dd','1111']
a(1)